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Abstract 

The dynamics of quantum expectation values is considered in a geometric setting. First, ex¬ 
pectation values of the canonical observables are shown to be equivariant momentum maps for the 
action of the Heisenberg group on quantum states. Then, the Hamiltonian structure of Ehren- 
fest’s theorem is shown to be Lie-Poisson for a semidirect-product Lie group, named the Ehrenfest 
group. The underlying Poisson structure produces classical and quantum mechanics as special 
limit cases. In addition, quantum dynamics is expressed in the frame of the expectation values, 
in which the latter undergo canonical Hamiltonian motion. In the case of Gaussian states, ex¬ 
pectation values dynamics couples to second-order moments, which also enjoy a momentum map 
structure. Eventually, Gaussian states are shown to possess a Lie-Poisson structure associated 
to another semidirect-product group, which is called the Jacobi group. This structure produces 
the energy-conserving variant of a class of Gaussian moment models previously appeared in the 
chemical physics literature. 


1 Introduction 


The expectation value dynamics for quantum canonical observables Z = iQ,P) has always attracted 
much attention, especially in connection to the correspondence principle. In this context, quantum 
effects are produced by noncommutativity of the canonical observables, thereby leading to Heisenberg 
uncertainty. In turn, uncertainty is associated to distributional effects that are typically encoded 
in Schrodinger’s probability density \'ip{x)\^. Then, one is led to look at statistical moments in or¬ 
der to obtain information about macroscopic quantities. Ehrenfest’s equations for expectation value 
dynamics read 

ih^^{Z) = {[Z,H]), (I) 

where H is the quantum Hamiltonian operator and the average is computed with respect to the 
density matrix (that is, p = for pure states). Ehrenfest dynamics has recently been considered 
from a geometric perspective [16] and the present paper uses similar geometric methods to identify 
its Hamiltonian structure. 

The relation (1) can also be expressed in terms of the Wigner phase-space function W{q,p) 
associated to p [48]. This is conveniently written in Dirac notation, so that in one dimension the 
ket |x) is the eigenvector given by Q\x) = x\x), while the bra vector is dehned as {x\ := jx)^. Upon 
generalizing to three dimensions, the Wigner function is associated to p by the relations [50] 
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The mapping that gives the Wigner function from the density matrix (first relation above) is known 
as Wigner transform of p, while the inverse mapping (second relation above) is called inverse Wigner 
transform, or more often Weyl transform of W{q,p). These mappings generate a one-to-one corre¬ 
spondence between operators and phase-space functions. The Wigner function is a quasi-probability 
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distribution (i.e., it can be negative) and is sometimes simply called ‘Wigner distribution’. Its evolu¬ 
tion reads dtW = {{H,Wy\, where H{q,p) is the Wigner transform of H (also known as the Weyl 
symbol) and {{•, •}} = {•, is the Moyal bracket [33, 50] carrying the quantum noncommuta- 

tive deviation (of order 0{h?)) from the classical Poisson bracket {f,g}^ = dxf ■ dpQ — d^g ■ dpf. The 
explicit definition of the Moyal bracket operator is rather involved and is omitted here. However, we 
shall recall that, although it does not identify a Poisson structure (as the Leibniz product rule is not 
satisfied), the Moyal bracket defines a Lie bracket structure. In this framework, Ehrenfest’s theorem 
for the expectation values (C) = /Cl^(CT)d®C is taken into the form d{C)/dt = ({{C,H}}), where 
C = (Qj p) is the phase-space coordinate. At this stage, we recall another fundamental property of the 
Moyal bracket of two phase-space functions: whenever one of these two is a second-degree polynomial, 
the Moyal bracket drops to the classical Poisson bracket. Consequently, {{C,H}} = and the 

Wigner-Moyal formulation of Ehrenfest’s theorem reads 

which coincides with classical expectation value dynamics (as it arises from the classical Liouville 
equation dtQ = {H,g}^), except for the fact that the averages are computed with respect to the 
quantum Wigner distribution W (whose evolution accounts for quantum noncommutativity). 

In the particular case when H (or its Weyl symbol H{Q) is a quadratic polynomial, the Wigner 
equation coincides with the classical Liouville equation and it is solved by a Gaussian distribution 
on phase-space. In this case, an initial Gaussian evolves in time by changing its mean and variance, 
so that the latter moment tends to zero in the formal limit h ^ 0 (i.e. the Gaussian tends to a 
delta function or point particle). See [31] for a geometric treatment of semiclassical Gaussian state 
evolution. Moreover, for quadratic Hamiltonians the Ehrenfest theorem reduces to classical particle 
dynamics, that is 

z = JVzH(z), 

where J is the canonical Poisson tensor = {C*;C'^}c we have defined z := {Q). Then, one is led 
to conclude that the expectation values follow canonical Hamiltonian trajectories as long as an initial 
wavepacket keeps narrow in time, so that the expansion 

H{C) ^ H(z) + (C - z) • VH(z) + i(C - z) ■ VVH(z)(C - z) (2) 

is justified and ({(^,H}^) ~ JVH(z). However, this statement leads to interpretative questions [5] 
and there is no reason why the approximation above should hold in the general case. 

Indeed, when H is not quadratic, an initial Gaussian Wigner function will evolve into different (non 
Gaussian) phase-space profiles, due to the effect of the noncommutative terms in the Moyal bracket. 
The time by which these quantum corrections become important (so that the initial Gaussian profile 
spreads over macroscopic scales) is known as Ehrenfest time and is related to the possible breakdown of 
classical-quantum correspondence in quantum chaos [34]. Then, a useful way to quantify the classical- 
quantum differences due to noncommutativity is through the difference between solutions of the 
classical Liouville and quantum Wigner equations (both initiated in the same state), respectively. In 
recent years, this has been attempted by comparing the dynamics of quantum and classical statistical 
moments [4, 12], that is the moments of the quantum Wigner function and the moments of the 
classical Liouville distribution, respectively. Still, the relation between expectation value dynamics 
and canonical Hamiltonian motion remains unclear and the present paper aims to shed some light on 
this point. For example, in [10] (see equations (6.2) therein) the expectation values z = {0 are shown 
to obey canonical Hamiltonian dynamics (even for non-Gaussian states) when the energy is expressed 
in terms of moments ((^ — z)^). However, this result still lacks a more fundamental description in 
terms of quantum state dynamics. 

Without entering further the difficult questions concerning the interpretation of Ehrenfest’s equa¬ 
tions, this paper unfolds their geometric properties and expresses the quantum evolution in the phase- 
space frame co-moving with the expectation values z(f). More particularly, if one writes the total 
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energy h{W) = J W{C,t) H{Q d®C terms of the Wigner function W (C, t) and the expectation values 
z (that is h{W) = h{W, z)), introducing the fluctuation variables C = C — z so that W (C, t) = W (C, t), 
yields the Ehrenfest equations in canonical form (see Section 3) 

z = JVzh(lE, z), (3) 

where h{W, z) = h{W, z) coincides with the total energy of the quantum system. As a new result, 
this paper also identifies the dynamics of the quantum state W{C,t) accompanying the evolution of z. 
Interestingly enough, equation (3) is identical in form to the equation of motion of a classical particle in 
the Ehrenfest mean-field model of mixed classical-quantum dynamics [19, 49]. To avoid confusion, it is 
important to point out that (3) does not mean that the expectation values follow classical trajectories, 
as they would be obtained by the classical limit of a quantum system. Indeed, while the canonical 
structure implies Hamiltonian trajectories of classical type, these trajectories do not coincMe with 
those of the classical physical system, which in turn would be obtained upon replacing h{W, z) by 
the classical Hamiltonian function h(z). This fact avoids the possibility of contradiction in equation 
(3), which is to be coupled to the evolution of the relative quantum distribution W. We address the 
reader to [15] for further remarks on the geometric analogies and differences between classical and 
quantum dynamics. 

In this paper, we shall focus on the geometric nature of the classical-quantum coupling that emerges 
from expectation value dynamics. The first goal is to present (in Section 2) a new formulation 
of Ehrenfest’s theorem in terms of a classical-quantum Poisson bracket that couples the classical 
canonical bracket to the Poisson structure underlying quantum dynamics. This classical-quantum 
bracket naturally incorporates classical and quantum mechanics as special cases. However, as it was 
briefly mentioned earlier, this work is not meant to provide a new interpretation of the classical- 
quantum correspondence. Rather, one of the targets of this paper is to express (in Section 3) the 
quantum dynamics in the phase-space frame co-moving with the expectation values, so that Ehrenfest 
equations possess a standard canonical form independently of the quantum state (not necessarily 
Gaussian). This canonical structure was found in [10] in the context of quantum cosmology: in 
this case, equation (3) is accompanied by the evolution of the moments = {{(^ — z)*'). In the 
present paper, we replace the moment hierarchy by the evolution of the quantum state IT(C,t) in 
the co-moving frame, where the entire Ehrenfest theorem is rewritten explicitly. The second part 
of the paper specializes this construction to consider the evolution of Gaussian states, as they are 
identified with Gaussian Wigner functions on phase-space. As a result, a Poisson bracket structure is 
provided in Section 4 for the Gaussian moments (C) and (CO- These quantities enjoy an (equivariant) 
momentum map structure that confers them a Lie-Poisson bracket [25, 27, 32]. Then, this bracket 
is analyzed in detail and related to the role of metaplectic transformations for the particular case of 
Gaussian wavepackets. Eventually, this construction is applied to provide energy-conserving Gaussian 
closure models. More particularly, the characterization of Gaussian state dynamics enables providing 
energy-conserving variants of previous Gaussian moment models [39, 41] in which the role of energy 
conservation has previously posed some issues [42]. As a general comment, we wish to emphasize 
that this paper does not deal with the convergence issues that may emerge in quantum mechanics. 
In particular, probability densities are assumed to decay sufficiently fast so that expectation values 
converge at all times. 


2 Hamiltonian structure of Ehrenfest’s theorem 

We start this Section with a mathematical result that lies at the basis of equation (3): expectation 
values are momentum maps for the standard representation of the Heisenberg group ~ x 

on the space of wavefunctions. For completeness, we recall the multiplication rule for the Heisenberg 
group: 

(Zi, (pi){z2, ip2) = (zi+Z2,ipi+^l-^Zi- IZ2 
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A momentum map is, generally speaking, given by the generating function of the canonical transfor¬ 
mation associated with a certain symmetry group G. For a symplectic vector space {V,u}) carrying a 
(simplectic) G—representation, the momentum map J{z) is given by 2{J{z),^) := z), for all 

z & V and all ^ G g. Here, g denotes the Lie algebra of G, Cv{z) denotes the infinitesimal generator 
of the G—representation and (•, •) is the natural duality pairing on g* x g (here, g* is the dual vector 
space of g). While the reader is referred to [25, 27, 32] for further explanations, here we specialize 
this definition to the Hilbert space of wavefunctions, endowed with the standard symplectic form 
= 2hlm{'tpl\^p2) (here, {ipi\ip 2 ) ■= J'4’i{x)'ip2{x) d^x). Indeed, the symplectic structure of the 
Hilbert space of wavefunctions is of paramount importance for the results in this paper. A momentum 
map associated to a unitary G—representation on the quantum state space is the map J(V') G g* given 
by (J('0),^) := h , where is the infinitesimal generator of the group action (with 

g (z G). When G is the whole group of unitary transformations, the corresponding momentum map 
is given by J(V') = [15]. For the Heisenberg group of phase-space translations, the action is 

the usual displacement operator of coherent state theory [37]. Upon denoting z = (q,p), this action 
is given by 

{x) := - q), (4) 

so that the infinitesimal generator reads ^(V’) = —+ C ' and we have identified ^ = 

{C,4>) £ [l(M^"') (here, f}(M^"') ~ denotes the Heisenberg Lie algebra so that the pairing (•,•) 

on X t)(M^"') is the standard dot product on Then, upon computing (f^(V’)IV') = 

/i“^((/>||V'|P + C ■ we obtain the momentum map 


so that, up to applying J, the momentum map associated to the Heisenberg group representation 
identifies the expectation values of canonical observables. This momentum map is equivariant, that is 
Ad*^ (^) J(V') = J(4*(z,¥>)('*/')) (here. Ad* denotes the coadjoint representation for the Heisenberg group 
'H(R^”)), and thus it is a Poisson map (see Section 12.4 in [32]). Most importantly, this map takes 
the Poisson structure on the symplectic space of wavefunctions into the Lie-Poisson structure on 
f}(R^”')*, i.e. the canonical Poisson bracket {•,•}£ (up to a multiplicative normalization factor ||V’||^)- 
Notice that the momentum map J(V') induces an equivariant momentum map d{p) = (JI(/9|Z), Tr(/?)) 
on the space of density matrices(here, Tr denotes the standard trace of trace-class operators): in this 
case, the infinitesimal action of the Heisenberg group reads ^{p) = -^Z^p], For convenience 

purposes, the remainder of this paper identifies quantum states with either density matrices or their 
Wigner transforms (unless otherwise specified). As we shall see, this leads to major simplifications in 
the treatment of Gaussian quantum states. 

The momentum map property of quantum expectation values suggests looking for the Poisson 
bracket structure of Ehrenfest’s equations. We recall that the latter have to be accompanied by the 
evolution of the quantum state, this being given by a wavefunction, a density matrix or its Wigner 
function. In order to find the Poisson bracket for the expectation values, we start with the following 
Poisson bracket [43] for the quantum Liouville equation ihp = [H,p]: 




—ih ^ 


^1\ 

_5p'5p\/' 


( 5 ) 


where [■, •] denotes the standard commutator and / and g are formally defined as function(al)s on 
the space of Hermitian operators. This Poisson bracket returns the quantum Liouville equation as 
p = {p, d/i/dp}, with the Dirac Hamiltonian functional h{p) = {H) (total energy). Here, 6h/6p is a 
Hermitian operator and we identify vector spaces of linear operators with their dual spaces by using 
the pairing {A,B) = Re(A|H), where {A\B) = Tr(A^H) is the natural inner product. Since Ehrenfest 
dynamics advances both expectation values (denoted by z = (Z)) and the quantum state (identified 
with the density matrix p), we allow to consider functionals of the type f{z,p). In practice, this 
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means that in (5) we allow for / to depend on p both explicitly and through the expectation value 
{Z) = Tr(/9Z), so that f{p) = f{p,z). Then, one can use the chain rule to write 


-^fip) = ^ p) + 

op oz dp 

and replacing this in (5) yields 

5p'5p\l' 


{/, 9 }(/ 5 ,z) = {Lg}c + ih ^ 


7 

5z ’ 5p 


A 7 A 

6z ' 5p 


— ih ^ ( 


( 6 ) 

( 7 ) 


where we have used the canonical commutation relations. The term {f,g}c is multiplied by the 
normalization constant Tr(p) and in (7) we have set Tr(p) = 1, which is the standard normalization 
of the density matrix. Without this particular choice, the change of variables z i—>■ Jz and p !->■ —ihp 
takes the Poisson bracket (7) into the (minus) Lie-Poisson [25, 27, 32] on the semidirect-product Lie 
algebra (Du(^) (here, u{J^) is the Lie algebra of skew-Hermitian operators on J^), endowed 

with the Lie bracket 




h(K2'i)(Du(jr) 


^0, —• JI^25 ^[C 2 • 


-1^,6]+ [6,6]). 


Here, we recall the general form of a (±)Lie-Poisson bracket, i.e. {f,g}±{p) = ±{p,[6f/Sp,6g/Sp]) 
(where // € g* is an element of some dual Lie algebra): the Lie-Poisson bracket (7) is obtained 
by choosing the minus sign and by specializing to the case g = This semidirect- 

product Lie algebra has been recently used in the variational formulation of the Ehrenfest theorem 
[11], while it is shown here to emerge naturally from the momentum map structure of expectation 
values. Indeed, the occurrence of the canonical Poisson bracket is due to the fact that expectation 
values are Poisson momentum maps that take (5) into the classical canonical structure. We shall call 
the Lie algebra the Ehrenfest algebra and its underlying Lie group 77(M^’^)(s)7/(J^) 

the Ehrenfest group (here, U{J^) denotes the group of unitary operators on Jif). The construction 
of the Ehrenfest group uses the group homomorphism provided by the representation (4), as it was 
presented in detail in Section IV.B of [11]. The bracket (7) will be called Ehrenfest bracket: this is the 
first example of a classical-quantum bracket that couples the canonical Poisson bracket underlying 
classical motion to the Lie-Poisson bracket (5) underlying quantum Liouville dynamics. However, 
notice that this Poisson bracket does not model the correlation effects occurring in the interaction of 
quantum and classical particles. Indeed, Poisson bracket structures modeling the backreaction of a 
quantum particle on a classical particle have been sought for decades and are still unknown despite 
several efforts [2, 3, 13, 29, 38, 40]. Rather, the Ehrenfest bracket governs the classical-quantum 
coupling (middle term in (7)) between expectation value dynamics (first classical term in (7)) and 
quantum state evolution (last term in (7)) for the same physical system. 

At this point, in order to write down explicit equations of motion, one has to find the expression 
of the total energy h{p) = {H) in the form h{p) = h{p,z). More particularly, we ask that h{p,z) 
is linear in p, so that it can still be expressed as an expectation value. For example, the kinetic 
energy in one spatial dimension can be rewritten by using the relation (R^) = + {{P — p)^), with 

p = (P). In more generality, one can Taylor expand the original Hamiltonian operator (or simply 
parts of it) around the expectation values. This leads to an expression of the total energy (H) of the 
form h{z,p) = {Hcq{z)), where Hcq{z) = 5h{z, p)/5p is a classical-quantum Hamiltonian operator 
that depends on the expectation values z (and not explicitly on p). For example, in the simplest case 
of a free particle in one spatial dimension one writes Hcq{z) = (p^/2)l -|- (P — pl)^/2. Evidently, 
analogous expressions hold for linear systems, such as the harmonic oscillator. Notice that Hcq{z) 
differs from the original quantum Hamiltonian 77, although they generate (by definition) the same 
total energy (77 cq(z)) = (77). Then, upon using the total energy expression h = (Hcq), the Poisson 
bracket (7) produces the equations 

z = JVz(Pcq) - ih-\[Z, Hcq] ) , 

ihp = [77 cq + Vz{Hcq) ■ Z, p]. 
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( 8 ) 

(9) 










Here, it is important to emphasize that the relation z = can be used only after evaluating the 
derivatives Vz(-Hcq)) which are to be computed by keeping z and p as independent variables (for 
example, one has Vp[p{P)) = (P) = p). 

Equations (8) and (9) were obtained in [11], upon postulating a specific variational principle, 
based on analogies with the variational structure of the Ehrenfest mean-field model of mixed classical- 
quantum dynamics. Then, the action principle postulated in [11] is justified here in terms of its 
corresponding Hamiltonian structure. More importantly, we have shown how these equations are 
totally equivalent to the quantum Liouville equation. Indeed, these equations were derived from the 
Poisson bracket (5) for the density matrix evolution, without any sort of assumption or approximation: 
the only step involved was rewriting the total energy as h{p) = h{p, z), which is no loss of generality as 
long as one term in the original Hamiltonian H can be expanded around the expectation values. Notice 
that, since equation (9) is indeed consistent with the ordinary Liouville equation, the classical-quantum 
Hamiltonian is constructed in such a way that it may differ from the original quantum Hamiltonian 
H only by a phase factor, so that H = Hcq + '^z{Hcq) ■ Z + 4>1 (upon dropping the phase factor, 
this is simply a consequence of the chain rule relation (6)): since by construction {Hcq) = {H), then 
(j) = —z • ^z{Hcq) and the two Hamiltonians are related by 

H = Hcq + {Z-z)- V^{Hcq) . 

Upon fixing a certain quantum system with Hamiltonian this is a consistency relation that is 
satisfied by the classical-quantum Hamiltonian. On the other hand, fixing a specific classical-quantum 
Hamiltonian, the above relation gives the corresponding quantum Hamiltonian operator. 

Depending on the specific form of the classical-quantum Hamiltonian, equations (8) and (9) allow 
for two different limits. Indeed, while Ehrenfest’s theorem is obtained in the case '^z{Hcq) = 0 (that 
is, in the purely quantum case when the Hamiltonian operator is not written in terms of expectation 
values and Hqq = H), the phase-type operator Hcq = h(z)l yields another limit in which expectation 
values follow classical particle trajectories, i.e. z = JVh(z). As noted in [11], the latter case takes the 
quantum equation (9) into 

iHp = [Vzh ■ Z,p] , 

which determines how the quantum state (together with Heisenberg’s uncertainty) is carried along 
the classical trajectories. The above Liouville equation is associated to a quantum Hamiltonian 
H = h(z) + {Z — z) ■ Vh(z) that is linear in Z and as such it generates coherent state evolution 
[37] of an initial pure state po = V’oV’o- Notice that Wigner-transforming the above quantum equation 
of motion yields the classical Liouville dynamics 

dtWiC) = {V,h-C,W{C)}c, (10) 

so that noncommutative quantum effects are absent, as expected in coherent state dynamics. On 
one hand, this comes as no surprise, since coherent states are widely known to be classical states in 
quantum optics. (In quantum optics, the word ‘classical’ does not refer to the motion of a physical 
classical particle, which would require a delta-like solution of (10) and is recovered only in the formal 
limit ^ ^ 0 of Gaussian solutions of the type (16)). On the other hand, this yields a suggestive 
picture in which classical trajectories carry a coherent quantum state, whose dynamics decouples 
from classical motion. 

While equations (8) and (9) recover Ehrenfest’s theorem as a special case, it is important to remark 
that they carry a redundancy, in the sense that equation (8) is simply the expectation value equation 
associated to (9) (if the latter is interpreted as a nonlinear nonlocal equation). This redundancy is not 
new, since it already occurs in Ehrenfest’s original equations. However, the redundancy of equation 
(8) can be eliminated by expressing the quantum dynamics in the frame of the expectation values. 
As we shall see, this operation separates the expectation values from the fluctuations arising from 
quantum uncertainty. While doing this is not generally trivial when using wavefunctions or density 
matrices, it becomes rather straightforward when using Wigner’s phase-space description. We remark 
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that changes of frames in the configuration space for quantum dynamics were studied in the past 
[30, 44], while changing the phase-space frame comes here as a new concept. This is the topic of the 
following Section. 


3 Quantum dynamics in the frame of expectation values 


As anticipated above, the splitting between quantum averages and fluctuations can be conveniently 
performed by writing the quantum dynamics in the frame of the expectation values. For this purpose, 
it is convenient to work in the Wigner phase-space formalism. Then, the equations (8) and (9) become 

dW 

z = JV.(//cq) + ({C, Hcq}) , = {{Hcq, VF}} + V,(i7cQ) • {C, , 

where Hcq{z,C) is the Weyl symbol of the classical-quantum operator Hcq{z). At this point, we 
perform the change of variables C = C + z, so that the Weyl symbol of the Hamiltonian is rewritten 
as Hcci{z,C) = H{z,C) and analogously W{C,t) = W{C,t). No confusion should arise here from 
the use of the tilde notation, which differs from that used in the previous Section. Notice that this 
change of variables (or change of frame) preserves the canonical commutation relations and thus it 
does not affect Heisenberg’s uncertainty principle. This change of variables returns (3), which then 
holds regardless of the specific type of quantum state. Also, the Wigner equation above transforms 
into dtW = {{H, W}} — (V^H) ■ {(^, W}. Notice, the same result can be obtained in terms of Poisson 
structures by changing the variables in the Ehrenfest Poisson bracket (7), after expressing the latter 
in terms of Wigner functions. This can be done by considering the Lie-Poisson bracket for the Wigner 
equation [8] (recall the notation (A(C)) = /kF(C)A((^) dC) 


and by following similar steps as in the previous Section to obtain the Ehrenfest bracket (7) in the 
form 


{f,g} = {/,fi'}z + 


^ A 

l(5z ■ 


^ ^ ILX 

5z 5 wL 


+ 


IIL AT 

XSW 5WI 


( 12 ) 


Here, we have denoted {f,g}z = Vz/ ■ IVzg and analogously for while the Moyal bracket 

operates only on the phase-space quantum coordinates (so that {{/,<?}} = {{/jS}}^)- Notice the abuse 
of notation, as / and g are very different in (11) and in (12) (compare with (5) and (7)). Then, 
changing variables to C = ^ — z, so that /(IF, z) = f(lF,z), yields the relations 


1L = ^ //il = A\\ 

dW 5vf’ [[SW'6W }} IUif’(51f//’ 

so that the Ehrenfest bracket becomes 


6z 6z 


IF Vj dC, 


{f,g}(lF,z) ={f,g}z + /|c,^| 

\ I oVF J ^ 



6f_ 

6W' SW 


(13) 


where we used and the notation (F(C)) ;= / 1F(C) F(C) d®/. Here, a useful feature 

of the Moyal bracket is the permutation property f a(C){{t>(C), c(C)}} d®C = f c(C){{«(C)) &(C)}} for 
any three functions a,b,c G C'°°(M®). Eventually, the total energy expression h(z, W) = {Hcq{z, C)) = 
{H{z,(^)) = h(z,lF) yields the coupled equations 


z = {z, {H)}z , = {{H, IF}} + {C • J({C, H}^), TF}^-. (14) 
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Here, (^) = 0 is verified to be preserved in time as the zero-level set of the expectation value momentum 
map (see [46] for an analogue of this in classical moment dynamics). Then, the redundancy occurring 
in (8) and (9) has now been eliminated, since the expectation value equation is no longer computed 
from the quantum Wigner-Moyal equation, which actually encodes the quantum deviations from the 
mean variables z. On the other hand, as already found in [10], the first equation in (14) implies 
that expectation values evolve along canonical Hamiltonian trajectories with Hamiltonian (H). This 
expectation value dynamics is different from classical Hamiltonian motion (one of the limits treated 
in the previous Section), as it carries higher moments of the type (CC • • • C) (upon assuming that H is 
analytic) that are responsible for quantum deviations from the classical physical trajectory [10]. The 
latter still emerges when H(z,(^) is linear in the deviation coordinate which is the case of coherent 
states. In this case, the relative distribution does not evolve in time and the relation dtW = 0 is 
simply equation (10) rewritten in the frame of expectation values. 

Once quantum dynamics is expressed in the frame of expectation values, one can continue to com¬ 
pute moments • • • C) of the deviation coordinate. The identification of the Hamiltionian structure 
for these moments was carried out in [10], in the context of quantum cosmology [9]. This Hamiltonian 
structure was recently compared [12] to that of classical moments of the Liouville equation (see also 
[26] for the classical case). The main difference between the results in [10] and those in the present 
paper is that here the moment hierarchy is replaced by the explicit evolution of the quantum state 
in the phase-space frame moving with z{t). For example, in terms of the density matrix p (as it is 
obtained by the Weyl transform of lF(C,t)), equations (14) are replaced by 

i = {z,{Hcg)h, ih^=[HaQ,p] + 'j:{[HcQ,Z])-[SZ,p] (15) 

where Z = Z — z and Hcq{z) is the operator obtained as the Weyl transform of H{z, (with respect 
to C)- This description has the advantage of overcoming the moment truncation problem, while it 
requires finding solutions of a nonlocal nonlinear quantum equation. 

Notice that, unlike the moment approach, the present method does not require the Weyl symbol 
of the original Hamiltonian to be analytic. As an illustrative example in one dimension, here we 
consider a unit mass subject to a step potential, so that the phase-space Hamiltonian is H{(^) = 
p'^I2 + pQ{q), where ^ is a physical parameter, 0 denotes Heaviside’s step function, and we have 
used the notation ^ = (q, p). Upon introducing the deviation coordinate C = C ~ 2(^)) writes 
H{z -I- (^) = {p + p)^/2 -b p0(g -b q) with the notation z{t) = {q{t),p{t)). Then, the total energy 
h = {H{z + ^)} reads 

h(z,W) = ^(p^ + (p^})+p(e(q+q)} 

= ^ // + 

and thus S^jSW = (p^ -b p‘^)l2 -b pQ{q + ^ = H{z,^) is the effective Hamiltonian simh that h = 
(Ff(z-bC)) = {H). In this case, the function H{z, (^) is not analytic and the total energy (H) cannot be 
expressed exactly in terms of the hierarchy of moments {CC • • • 0- However, in the present approach, 
one replaces moments by the explicit quantum evolution in the frame moving with the expectation 
values z(t). In this framework, the latter evolve according to the first equation in (14), that is 

q=p, p= -p jw {-q, p) dp. 

It is not surprising that the expectation values evolve under the time-dependent effective potential 
P fdpW{q, p), so that the step potential is modulated to allow for quantum tunneling. Notice 

that the term ff p^/2 W{q, p) dp dp in the total energy {H) does not contribute to expectation value 
dynamics. In addition, the quantum state evolution is given by the Wigner equation (14) for W{q, p, t). 


This is given by the following integro-differential equation 


dW 


+ fie{q + p) ,W }} = -p 


dW 

dp 


Jwi 


-q,p) dp, 


where a new nonlocal term emerging from the change of frame appears on the right hand side. If 
this term is set to zero, one obtains the ordinary quantum dynamics in phase space for a particle 
subject to a potential step in —q. While explicit solutions can be constructed by taking the Wigner 
transform of the usual wavefunction of a unit mass subject to a step potential, the detailed study of 
the nonlinear nonlocal evolution of W is note among the purposes of this paper. 

Notice, the whole treatment proceeds analogously for classical Liouville dynamics upon replacing 
Moyal brackets with Poisson brackets. Once more, this means that the essential difference between 
classical and quantum statistical effects lies in the noncommutative terms (higher order in h) of the 
Moyal bracket. It is useful to emphasize that the conserved total energy {H) does not depend explicitly 
on time, as it is readily seen by expanding 0 = d{H)/dt = dt{H) + z • Vz{H) and using the first of 
(14), so that z • 'S/z{H) = 0. Notice that the relation dt{H) = 0 can also be verified explicitly by 
computing dt{H) = dtW{Q) d®(^ = 0. Moreover, another consequence of the canonical nature 

of expectation value dynamics is the invariance of Poincare’s relative integral invariant: Jz • dz = 

§^{P) ■ d{Q) = const., for any loop 7 in phase-space moving with the Hamiltonian vector field 

{z,(^)}z. 

In the next Section, we focus on Gaussian quantum states, thereby restricting to consider only 
second-order moments. More particularly, we shall characterize the Hamiltonian structure of Gaussian 
state dynamics in terms of expectation values z and covariance matrix (CO- 


4 Hamiltonian structure of Gaussian quantum states 

Once the dynamics of Ehrenfest expectation values has been completely characterized in terms of 
Poisson brackets, one may consider higher order moments. The moment hierarchy does not close in the 
general case, although it is well known that it does for quadratic Hamiltonians. In the latter case, the 
moment algebra acquires an interesting structure, which is the subject of the present Section. Before 
entering this matter, we emphasize that quadratic Hamiltonians restrict to consider linear oscillator 
motion and so they are uninteresting for practical purposes. An interesting situation, however, occurs 
when the total energy {H) is computed with respect to a Gaussian state, so that higher moments 
are expressed in terms of the first two. This is the Gaussian moment closure for nonlinear quantum 
Hamiltonians. 

Gaussian quantum states (a.k.a. squeezed states in quantum optics) have been widely studied over 
the decades in many different contexts, mostly quantum optics and physical chemistry. Recently, ap¬ 
plications of Gaussian states in quantum information have also been proposed (see e.g. [1]). Generally 
speaking, a Gaussian quantum state is identified with a Gaussian Wigner function of the form 

^ vdrts(t) 

where is a normalising factor, z = {(^) is the mean and AS = h{{CC) ~ (0(0) is the covariance 
matrix. The question of which covariance matrices are associated to genuine Gaussian quantum 
states was addressed in [31] for wavepackets and in [22, 45] for more general mixed states. For 
example, when the Gaussian is too narrow, then the corresponding density matrix is not positive- 
definite and the uncertainty principle is violated. However, when the initial covariance matrix satisfies 
Heisenberg’s principle (so that (16) identifies a quantum state), the latter holds at all times as the 
quantum uncertainty is preserved by the Wigner-Moyal equation dtG = {{H, G}} for Gaussian states. 
We emphasize that the expression (16) incorporates the Wigner transform of Gaussian wavepackets 
[23, 24] as a special case. 
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If the linear form w and the quadratic form S defining the quadratic Hamiltonian H = (^-SC+'w-(^ 
are functions (possibly nonlinear) of the first and second order moments (C) and (CC); then a Gaussian 
initial state will remain a Gaussian under time evolution by changing its mean and variance. In more 
generality, a nonlinear (analytic) Hamiltonian will produce a total energy (H) = f G(^}H(C} dC that 
can be expressed entirely in terms of first and second order moments, i.e. 


(H) = h((0,(CC)). 


At this point, it is useful to restrict the Ehrenfest bracket (7) (or its phase space variant (12)) 
to functionals of Gaussian Wigner functions G(^), depending only on the first two moments. For 
convenience, we shall denote X = {(^0/2, so that h{z,G) = H(z, A). The corresponding Poisson 
structure is easily found by using the chain rule relation 


6f /<5F 1 SF \ 

JG 2hAV ’ 

in (11). As a result, one finds 

{F. G)(z. X) = {F. G}. + z . ) + iv(x 

along with the following equations of motion for an arbitrary total energy H(z, A): 


SX^SX ~ 6X^6X I ’ 


5H 

z = {z, H}c + JI—z. 


A = 


6X 5X 


+ 




(17) 


(18) 


(19) 

( 20 ) 


The Poisson bracket (18) has appeared earlier in the literature [20, 28] in the context of classical 
Liouville (Vlasov) equations. This is no surprise, as second order moments are associated to quadratic 
phase-space polynomials, which do not involve the higher-order noncommutative terms in the Moyal 
bracket (Gaussian quantum states undergo classical Liouville-type evolution). Under the change of 
variables A AJ and z Jz, the above moment bracket was shown [20] to be Lie-Poisson on the 
Jacobi group 

= Sp(m2^)(D77(m2^) , 

i.e. the semidirect product of the symplectic group with the Heisenberg group. Here, the semidirect- 
product structure is defined by the following action of the symplectic group on the Heisenberg group: 
<l>5(z,(/3) = {Sz,ip), where S € Sp(M^”). Therefore, we conclude that the moments of any Gaussian 
quantum state evolve on coadjoint orbits of the Jacobi group. In particular, this means that the sym¬ 
plectic forms recently found to underlie the dynamics of Gaussian wavepackets [35, 36] are symplectic 
forms on coadjoint orbits of Jac(M^”'), which are determined by the usual Gasimir invariants [18, 26] 


C,(z,A) = lTr 


A--zz)J 


2 i 


where j = 1,2,3 and Tr denotes the matrix trace. For example, setting j = 1 and expanding yields 
(up to multiplicative factors) Ci(z, A) = det(2A — zz). 

These results come as no surprise. The relation between the Jacobi group and Gaussian states has 
been known for decades in the theory of coherent states [6, 7], under the statement that wavepackets 
evolve under the action of the semidirect product Mp(R^"')(D77(M^"') [17, 31], where Mp(M^”) is 
the metaplectic group (that is, the double covering of Sp(M^"’)). Since the Lie algebra mp(M^”') 
of the metapletic group is isomorphic to that of the symplectic group (denoted by sp(]R^")), then 
joc(M^") ~ mp(M^"') (D 1)(M^"'). In the present treatment, the symplectic group replaces the metaplectic 
transformations because we have identified quantum states with Wigner functions (which can account 
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for pure as well as mixed states) rather than wavefunctions. Indeed, while the symplectic group does 
not possess a representation on wavefunctions, it does possess a natural action on the space of Wigner 
(phase-space) functions as in classical mechanics and one may avoid dealing with the metaplectic 
representation [31]. This leads to an action of the Jacobi group on the space of Wigner functions, 
which is given by 

{<^is,.,^){W)m = W{SC + z), (21) 

where 5 is a symplectic matrix and {z,(p) is an element of the Heisenberg group. This action is given 
by the pullback of the Wigner function 1H(C) by the phase-space transformation ^ i—>■ + 2 :. 

The emergence of a Lie-Poisson bracket for the first and second moments is also not surprising. 
Indeed, this is due to the fact that the moment triple {{l),z,X) is itself an equivariant momentum 
map for the action (21); see Section III.C of [20]. (Here, we have formally denoted (1) = J’W(C)dC). 
Then, Gaussian Wigner functions are identified with the moment couple {z,X) and this identification 
enables the description of Gaussian state dynamics in terms of coadjoint orbits. 


5 Gaussian moment models and energy conservation 


Notice that one can rewrite the above dynamics in terms of the covariance matrix S = 2X — zz. This 
is easily done by restricting the bracket (13) to functions of the type H(z, S). This process yields the 
direct sum bracket 


{f,9} ={/,5}z+2Tr(s 








( 22 ) 


This Poisson bracket produces the following equations for a total energy of the form h = h{z,T,): 


z = {z,h}c, 




(23) 


For example, in the particular case when the total energy {H) = JG{C)H{C) d( is approximated 
by using (2), these equations recover the dynamics (12) and (13) in [21], suitably specialized to 
Hermitian quantum mechanics. The explicit comparison of the equations (23) with those obtained 
by Heller for Gaussian wavepackets (see [23, 24] and subsequent papers by Heller on the same topic) 
requires expressing the covariance matrix as a function on the Siegel upper half-space [36, 35]. The 
corresponding geometric description is the subject of ongoing work. 

Equations (23) can be directly applied to modify certain moment models that have previously 
appeared in the chemical physics literature [39, 41, 42]. This class of models suffers from lack of 
energy conservation in the general case [42] , with possible consequent drawbacks on the time evolution 
properties. In references [39, 41, 42] and related papers on the same topic, a class of moment models 
was developed by adopting a Gaussian moment closure on the equations of motion for the expectation 
values (Z) and {ZZ). More particularly, Gaussian closures of the type {ZZZ) ~ ‘i{ZZ){Z) — 2{Z)‘^ 
(and similarly for higher order moments) are used in the equations 

ihj^{Z) = {[Z,H]), {ZZ) = {[ZZ,H]). 

The moment closures are performed after replacing the Hamiltonian operator H by its Taylor ex¬ 
pansion around the expectation value z = {Z) (similarly to (2)). When this Taylor expansion is 
truncated to third order, one recovers the Gaussian moment closure previously introduced in [47]. 
However, in the general case of a higher-order expansion, adopting the closure directly in the Gaus¬ 
sian moment equations may break energy conservation [42]. Indeed, while odd-order expansions of H 
preserve energy conservation, this does not hold for truncations of even order greater than 2. More 
particularly, as noticed in [42], in order to conserve the total energy, the Taylor expansion used in 
the equations should be terminated after an odd derivative. The framework presented in this Section 
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provides a solution to this problem. Indeed, once the closure has been performed in the expression of 
the total energy H = (H), the Hamiltonian moment equations (19)-(20) are uniquely determined. As 
a consequence of the Poisson bracket structure (18), these equations generally differ from those in [42] 
by exactly the odd derivatives that ensure conservation of both the total energy and the determinant 
of the covariance matrix. More particularly, if one expands the Hamiltonian H = P^/2 + V{Q) up 
to fourth order around the expectation values {Z), the second parenthesis in equation (20) produces 
conservative terms consisting of the fifth-order derivatives in the following equations of motion (in 
standard expectation value notation): 

|(C)=(P>. 

A(p) = - H(1)((Q)) - ^H(^)((Q)) ((Q2) _ (Q)2) _ 1 h(5)((Q)) ((Q2) _ (Q)2)2 ^ 
^(Q2)=2 (qP),, 

^(p2) = _ 2 H(i)((Q)) (P) - 2 H(2 )((q)) _ (g)(p)) _ H(3)((g)) (p)((g2^ _ 

- h(")((Q)) {{QP)s - {Q){P)){{Q^) - {Qf) -{P){&) - {Q?f , 

^(QP), = (P 2 ) - h(i\(Q)) (Q) - ((Q2 ) _ (q)2 ) _ iH(3)((g)) (g)((Q2) _ ^g)2)^ 

- - {Q?f - \y^^\{Q))m{Q^) - {Qf ?, 

where {QP)s denotes the expectation of the symmetrized product {QP + QP)j2. Dropping the 
fifth-order derivatives returns the non-conservative equations (15)-(19) in [42]. Therefore, for exact 
polynomial potentials V(Q) of fourth degree, energy conservation is not an issue because V^^\{Q)) = 0 
and the two models coincide. However, more general cases such as Morse-type potentials require extra 
care in dealing with the Gaussian closure. The implications of the energy-conserving terms in concrete 
physical problems will be the subject of future work. 


6 Conclusions and perspectives 

Based on the Hamiltonian Poisson bracket approach, this paper has unfolded the geometric properties 
of Ehrenfest’s expectation value dynamics. More particularly, the search for the Hamiltonian struc¬ 
ture of Ehrenfest’s theorem has produced a new classical-quantum Poisson structure that incorporates 
classical and quantum dynamics as special cases. The corresponding equations are Lie-Poisson for 
the Ehrenfest group This result was achieved by exploiting the momentum map 

property underlying Ehrenfest expectation values. Later, the Ehrenfest theorem was rewritten upon 
expressing quantum dynamics in the phase-space frame co-moving with the expectation values. This 
result extends previous work in quantum cosmology [9, 10] by avoiding moment truncation prob¬ 
lems. In the last part, the Poisson structures underling Ehrenfest theorem were restricted to consider 
Gaussian moment dynamics, as it arises from Gaussian Wigner functions on phase-space. Again, 
Gaussian moments enjoy a momentum map structure that confers them a Lie-Poisson bracket for the 
Jacobi group Sp(]R^”)(S)?^(]R^"'). This bracket structure enables providing energy-conserving variants 
of previous Gaussian moment models [39, 41] that were generally lacking conservation of energy [42]. 

This paper has shown that the use of Wigner functions and the properties of the Moyal bracket 
are particularly advantageous for studying expectation values. Then, combining Poisson brackets 
with momentum map structures unfolds the geometry underlying quantum dynamics. Eor example, 
momentum map structures may also appear in quantum hydrodynamics, where local averages (e.g. 
JW{q,p) d^p and fpW{q,p)d^p) are considered. In addition, although the present treatment did 
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not consider spin effects, these can be retained by including the expectation {S) of the spin operator 
S in the treatment of the classical-quantum bracket (7). Another open question concerns the ladder 
operator formulation of these results. This can be particularly advantageous for the study of coherent 
squeezed states in quantum optics and is currently under development by the authors. Future research 
directions will also involve the Heisenberg picture and the corresponding variational framework. 


Acknowledgments. The authors are indebted to Tomoki Ohsawa, Alessandro Torrielli and the 
anonymous referees for providing extensive and valuable feedback on these results. In addition, the 
authors are grateful to Dorje Brody, Prangois Gay-Balmaz, Darryl Holm, David Meier, Juan-Pablo 
Ortega and Paul Skerritt for several discussions on this and related topics. This work was par¬ 
tially carried out at the Bernoulli Center of the Swiss Federal Institute of Technology in Lausanne: 
C.T. acknowledges hospitality during the program “Geometric Mechanics, Variational and Stochastic 
Methods”. Financial support by the Leverhulme Trust Research Project Grant 2014-112, the London 
Mathematical Society Grant No. 31320 (Applied Geometric Mechanics Network), and the EPSRC 
Grant No. EP/K503186/1 is also acknowledged. 


References 

[1] Adesso, G.; Ragy, S.; Lee, A.R. Continuous variable quantum information: Gaussian states and 
beyond. Open Syst. Inf. Dyn. 21 (2014), n. 1 & 2, 1440001 

[2] Aleksandrov, I.V. The statistical dynamics of a system consisting of a classical and a quantum 
subsystem. Z. Naturforsch. 36a (1981), 902-908 

[3] Anderson, A. Quantum backreaction on “classical” variables. Phys. Rev. Lett. 74 (1995), n. 5, 
621-625 

[4] Ballentine, L.E. Dynamics of classical-quantum differences for chaotic systems. Phys. Rev. A 
65 (2002), 062110 

[5] Ballentine, L.E.; Yang, Y.; Zibin, J.P. Inadequacy of Ehrenfest’s theorem to characterize the 
classical regime. Phys. Rev. A 50 (1994), n. 4, 2854-2859 

[6] Berceanu, S.A.A. Coherent states associated to the Jacobi Group: a variation on a theme by 
Erich-Kdhler, J. Geom. Symmetry Phys. 9 (2007), 1-8 

[7] Berceanu, S.A.A. The Jacobi group and the squeezed states: some comments. AIP Conf. Proc. 
1191 (2009), 21-29 

[8] Bialynicki-Birula, L; Morrison, P.J. Quantum mechanics as a generalization of Namhu dynamics 
to the Weyl-Wigner formalism. Phys. Lett. A 158 (1991), 453-457 

[9] Bojowald, M. Quantum Cosmology. Lect. Notes Phys. 835. Springer, 2011. 

[10] Bojowald, M.; Skirzewski, A. Effective equations of motion for quantum systems. Rev. Math. 
Phys. 18 (2006), n. 7, 713-745 

[11] Bonet Luz, E.; Tronci, C. Geometry and symmetry of quantum and classical-quantum variational 
principles. J. Math. Phys. 56 (2015), n. 8, 082104 

[12] Brizuela, D. Statistical moments for classical and quantum dynamics: Formalism and generalized 
uncertainty relations. Phys. Rev. D 90 (2014) 90, 085027 

[13] Caro, J.; Salcedo, L.L. Impediments to mixing classical and quantum dynamics. Phys. Rev. A 
60 (1999), n. 2, 842-852 


13 


[14] Chernoff, P.R.; Marsden, J.E. Some remarks on Hamiltonian systems and quantum mechanics. 
Univ. Western Ontario Ser. Philos. Sci. 6c (1976), 35-53 

[15] Chruscihski, D.; Marmo, G. Remarks on the GNS representation and the geometry of quantum 
states. Open Syst. Inf. Dyn. 16 (2009), 157-177 

[16] Clemente-Gallardo, J.; Marmo, G. The Ehrenfest picture and the geometry of quantum mechan¬ 
ics. II Nuovo Cimento G 36 (2013), n. 3, 35-52 

[17] De Gosson, M. Symplectic geometry and quantum mechanics. Birkhauser Verlag, Basel, 2006. 

[18] Dodonov, V.V.; Man’ko, O.V. Universal invariants in quantum mechanics and physics of optical 
and particle beams. J. Rus. Laser Res. 21 (2000), n. 5, 438-464 

[19] Ehrenfest, P. Bemerkung iiber die angendherte Giiltigkeit der klassischen Mechanik innerhalb 
der Quantenmechanik. Zeitschrift fiir Physik 45 (1927), n. 7-8, 455-457 

[20] Gay-Balmaz, F.; Tronci, G. Vlasov moment flows and geodesics on the Jacobi group. J. Math. 
Phys., 43 (2012), n. 12, 123502 

[21] Graefe, E.M.; Schubert, R. Wave-packet evolution in non-Hermitian quantum systems. Phys. 
Rev. A 83 (2011), n. 6, 060101(R) 

[22] Gracia-Bondia, J.M.; Varilly, J.C. Nonnegative mixed states in Weyl-Wigner-Moyal theory. 
Phys. Lett. A 128 (1988), n. 1-2, 20-24 

[23] Heller, E.J. Time-dependent approach to semiclassical dynamics. J. Ghem. Phys. 62 (1975), 
1544-55 

[24] Heller, E.J. Glassical S-matrix limit of wavepacket dynamics. J. Ghem. Phys. 65 (1976), 4979-89 

[25] Holm, D.D. Geometric Mechanics - Parts I &: H. Imperial Gollege Press. 2nd edition, 2011. 

[26] Holm, D.D.; Lysenko, W.P.; Scovel, J.C. Moment invariants for the Vlasov equation. J. Math. 
Phys. 31 (1990), n. 7, 1610-1615 

[27] Holm, D.D.; Schmah, T.; Stoica, C. Geometric Mechanics and Symmetry. Oxford University 
Press, 2009 

[28] Holm, D.D.; Tronci, C. Geodesic Vlasov equations and their integrable moment closures. J. 
Geom. Mech. 1 (2009), n. 2, 181-208 

[29] Kapral, R. Progress in the theory of mixed classical-quantum dynamics. Annu. Rev. Phys. Ghem. 
57 (2006), 129-157 

[30] Levy-Leblond, J.-M. Galilei group and nonrelativistic quantum mechanics. J. Math. Phys. 4 
(1963), 776 

[31] Littlejohn, R.G. The semiclassical evolution of wave packets. Phys. Rep. 138 (1986), n. 4-5, 
193-291 

[32] Marsden, J.E.; Ratiu, T.S. Introduction to Mechanics and Symmetry, Springer, 2nd edition, 
2002 

[33] Moyal, J.E. Quantum mechanics as a statistical theory. Proc. Cambridge Phil. Soc. 45 (1949), 
99-124 

[34] Karkuszewski, Z.P.; Zakrzewski, J.; Zurek, W.H. Breakdown of correspondence in chaotic sys¬ 
tems: Ehrenfest versus localization times. Phys. Rev. A 65 (2002), 042113 


14 


[35] Ohsawa, T. Symmetry and conservation laws in semiclassical wave packet dynamics. J. Math. 
Phys. 56 (2015), 032103 

[36] Ohsawa, T.; Leok, M. Symplectic semiclassical wave packet dynamics. J. Phys.A.: Math. Theor. 
46 (2013), n. 40, 405201 

[37] Perelomov, A. Generalized coherent states and their applications. Springer, 1986. 

[38] Prezhdo, O.V. A classical-quantum bracket that satisfies the Jacobi identity. J. Chem. Phys. 124 
(2006), 201104 

[39] Prezhdo, O.V. Quantized hamilton dynamics. Theor. Chem. Acc. 116 (2006), 206-218 

[40] Prezhdo, O.V.; Kisil, V.V. Mixing quantum and classical mechanics. Phys. Rev. A 56 (1997), n. 
1, 162-175 

[41] Prezhdo, O.V.; Pereverzev, Yu.V. Quantized Hamilton dynamics. J. Chem. Phys. 113 (2000), n. 
16, 6557-6565 

[42] Prezhdo, O.V.; Pereverzev, Yu.V. Quantized Hamilton dynamics for a general potential. J. 
Chem. Phys. 116 (2002), n. 11, 4450-4461 

[43] Requist, R. Hamiltonian formulation of nonequilibrium quantum dynamics: Geometric structure 
of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy. Phys. Rev. A 86 (2012), 022117 

[44] Schmutzer, E.; Plebanski, J. Quantum mechanics in non-inertial frames of reference, Fortschr. 
Physik, 25 (1977), N. 1-12, 37-82 

[45] Simon, R.; Sudarshan, E.C.G.; Mukunda, N. Gaussian Wigner distributions: a complete char¬ 
acterization. Phys. Lett. A 124 (1987), 223 

[46] Tronci, C. Hamiltonian approach to hybrid plasma models. J. Phys. A 43 (2010), n. 37, 375501 

[47] Turner, R.E.; Snider, R.F. A phase space moment method for classical and quantal dynamics. 
Can. J. Phys. 59 (1981), 457-470 

[48] Wigner, E.P. On the quantum correction for thermodynamic equilibrium. Phys. Rev. 40 (1932), 
749-759 

[49] Wyatt, R. Quantum dynamics with trajectories: Introduction to quantum hydrodynamics. 
Springer. 2005 

[50] Zachos, C.; Fairlie, D.; Curtright, T. Quantum mechanics in phase space: an overview with 
selected papers. World Scientific, 2005. 


15 


